Long hold times in a two-junction electron trap 
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The hold time r of a single-electron trap is shown to increase significantly due to suppression of 
environmentally assisted tunneling events. Using two rf-tight radiation shields instead of a single 
one, we demonstrate increase of r by a factor exceeding 10^, up to about 10 hours, for a trap with 
only two superconductor (S) - normal-metal (N) tunnel junctions and an on-chip resistor R (R-SNS 
structure). In the normal state, the improved shielding made it possible to observe r ~ 100 s, which 
is in reasonable agreement with the quantum-leakage-limited level expected for the two-electron 
cotunneling process. 



Creating a quantum current standard based on 
frequency-driven single-electron transport has been one 
of the major goals of metrology for about 25 years 
In the 1990's, the efforts were focused mainly on metallic 
tunnel junction devices. Relative uncertainty of about 
10^^ was reached by the seven-junction electron pump 
for small currents at the level of a few pA [H . Recent ad- 
vances have raised the hope to increase the current 
by factor 10-100 towards a more practical level, renewing 
the interest to the accuracy issues. 

Photon (aka environmentally) assisted tunneling 
(PAT) 0] has been found to limit the transport accu- 
racy at the level far above theoretical predictions [1, Q . 
Recently, it was shown that PAT rate T in the super- 
conductor (S) - normal metal (N) hybrid turnstile Q 
can be reduced down to time-resolved single events by 
using on-chip noise filteringelements: a high-ohmic re- 
sistor (R-SNS turnstile [13, El) or a capacitively coupled 
ground plane [7l. [l2l . [l3j . The electron trapping properties 
of the pumps are important for electron-counting capaci- 
tance standards [l3| , and also a way to study the leakage 
mechanisms. In the preceding experiment with the R- 
SNS turnstile-based electron trapping circuit, Ref. [llj, a 
reasonably low value of T = t~-^ ^ (20 s)~^ was achieved 
in a strong Coulomb blockade. Here, r is the average 
hold time of an electron in the trap. 

In this Letter, we report on the very same R-SNS trap 
samples tested in a cryogenic setup with a different type 
of noise filtering: double rf shielding of the sample stage. 
We demonstrate more than 1000-fold improvement of r 
up to about 10 h. This result is comparable with the 
hold times previously reported for much more complex 
arrays with 7 or more tunnel junctions [isl [l6j or with 
4 junctions and a high-ohmic resistor connected in series 
[17| . With the superconductivity suppressed by a mag- 
netic field, tunneling rates down to the level expected for 
two-electron cotunneling (l8l - l20| were achieved. 

The equivalent circuit of the trapping device is shown 
in Fig. 1(a). The left junction of an Al-AuPd-Al SNS 
turnstile is connected to an on-chip Cr resistor with 



R ^ 100 kr2, and the right junction is terminated by a 
trapping island whose charge state n2 is monitored by a 
capacitively coupled SNS single-electron transistor (SET) 
electrometer operating in the voltage-bias regime. For 
the details of sample fabrication see Ref. [ll| . 

Two nominally similar samples, Sa and Sb, were mea- 
sured in two cryostats, Ca and Cb, respectively. Both 
sample stages were equipped with two nested rf radia- 
tion shields. In setup Ca, the signal lines were filtered 
by the combination of Thermocoax'^'^and powder filters, 
and in Cb only by Thermocoax"'"'^. The radiation shields 
of Ca were of higher quality, and unless otherwise stated, 
the reported results are from Sa/Ca- During the prepa- 
ration of this report, we have applied similar nested rf 
shielding in studies of sub-gap behaviour of SNS tran- 
sistors 1131, and we also became aware of the work of 



' I antti .kemppiiien@mikes .fi I 



Barcnds et al. |2l| in which similar techniques were used 
to lower loss and decohcrencc of superconducting quan- 
tum circuits. 

At low temperatures, T < 100 mK, tunneling in 
the zero-biased trap is blocked by the energy barrier 
formed by the quasiparticle excitation energy Aai = 
260 /xeV = feflX 3.0 K and charging energies Eci ~ 
eV [2 X (2Ct + Cg)] « fcsx 7.6 K > Ec2 « eV2Ct,.ap « 
/cbX 3.3 K (estimated for sample Sa) of an electron lo- 
cated in the middle and in the trapping island, respec- 
tively. Here Ct is the tunnel capacitance assumed equal 
for both junctions, Cg is the gate capacitance of the turn- 
stile, and Ctrap is the total capacitance of the trapping 
island. The corresponding stability diagram is shown 
in Fig. 1(b) with the hold times r calculated using the 
orthodox theory of single-electron tunneling [l[ for four 
charge states of the trap in thermal equilibrium. 

Electrons can be loaded or unloaded from the trap by 
ramping the bias voltage Vb and thus lowering the en- 
ergy barrier for one of the two tunneling directions. Due 
to the energy barrier, the bias voltage dependence of the 
trapped charge n2 monitored by electrometer current lei 
appears hysteretic, see Fig. 1(c). The maximum of the 
barrier height and, thus, the width of the hysteresis loop 
can be controlled by the gate voltage Vg. Figure 1(c) 
shows several ramps of Vb encircling the widest loop 
which corresponds to an integer value of Ug = CgVg/e, 
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FIG. 1. (Color online) (a) Circuit diagram of the trapping de- 
vice. The number of electrons on the normal-metal (AuPd) 
middle island, ni, and on the superconducting (Al) trap is- 
land, n2, can be controlled by the voltages Vb and Vg. (b) 
Simulated contour plot of r at 120 mK for four charge con- 
figurations (ni,n2) as a function of (7ig,Vb), rig = CgVg/e. 
Bright and dark areas denote long and short values of r, re- 
spectively. In the bistable areas, the largest values of r are 
plotted. Inclined rectangles around each stable charge state 
correspond to r = 1 s escape thresholds. Two inclined solid 
lines show the middle points of the hysteresis loops for the 
states (0,0)-(0, 1) and (1,0)-(1,1). (c) Repeated Vb sweeps 
around the widest hysteresis loop of about 1.4 mV, empha- 
sized by the rectangle corresponding to the wide grey vertical 
line in (b). The red dashed line in the middle denotes the 
symmetric case, with equal values of r for both charge states 
involved. Inset: Electrometer signal Id as a function of gate 
charge rig^ei = Cg,ciV^,ci/e. The circle denotes a typical oper- 
ating point, (d) Random switchings of the trap between two 
charge states. The trap is tuned to exhibit experimentally 
convenient hold times by setting the width of the hysteresis 
loop to about 100 fjY. The central panel shows the symmetric 
case with Vh set in the middle of a hysteresis loop (AVb — 0), 
while the upper and lower traces are for AVb = t7 which 
favor one of the states. 



see Fig. 1(b). Due to finite T and external excitations 
particularly addressed in this work, the observed switch- 
ings between charge states include a random component, 
leading to an uncertainty of the loop width and to the 
two- level fluctuator behavior depicted in Fig. 1(d). 

To support the conckisions of the long hold times 
shown below, we studied r as a function of the energy 
barrier (electrostatic energy change) of tunneling, AE, 
which can be controlled by voltages Vb and Vg. Two-state 
switching traces were analyzed statistically for different 
values of AE. Figure 2(a) shows typical dependencies 
of T on the deviation of bias voltage from the middle 
of the hysteresis loop (AVb) for three different values of 
the electrometer current When AVb = 0, AE is the 
same for the two values of n2. A shift in AVb shifts 
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FIG. 2. (Color online) (a) Average hold time r of two charge 
states as a function of bias voltage deviation AVb from the 
middle of the hysteresis loop measured at three different levels 
of the electrometer current: I^i < 100 pA (black squares). 
Id ~ 300 pA (blue circles), and Id ~ 750 pA (red triangles). 
Open and filled symbols of the same type correspond to the 
two alternating charge states within the trace: n2 = and 
712 = 1, respectively [see Figs, l(c-d)]. For each Id, Vg was 
adjusted such that r « 1 s at AVb = 0. (b) Summarizing plot: 
hold time r vs. energy barrier AE. Green solid line shows the 
slope corresponding to thermal activation at T* — 120 mK. 
Symbols are as in (a). Black diamonds and pentagons were 
measured at Id < 100 pA with a wider hysteresis loop than 
black squares. 



proportionally both values of AE, but in opposite direc- 
tions. The observed linear dependency of log(r) on AE 
is typical for thermal activation (TA, cf, e.g., Ref. jT5|): 
T oc cxp{~ AE / {k-sT*)) , where the effective temperature 
T* characterizes the spectrum of fluctuations driving the 
state switchings. 

To study the crossover from TA to PAT, we present a 
plot in Fig. 2(b) summarizing the hold time data r ver- 
sus AE for different electrometer currents Ie\. For this 
approximative plot, we used the following conventions: 
(i) The origin of the AE = axis was chosen arbitrarily 
as the point where r « 0.1 s. (ii) The proportional- 
ity factor between AE and AVb is based on the fitted 
stability diagram. (Hi) For the curves with the smallest 
Id (yielding the longest hold times), we combined the 
piccewise-continuous data sets measured at different Vg 
and arranged them along the Ai?-axis in Fig. 2(b) to 
ensure their smooth alignment. This was done because 
collecting the statistics for a short-living state would take 
impractically long if following a large hysteresis loop: 
switching back from the long-living counter state would 
take hours. 

Although Fig. 2(b) is approximative, it allows us to 
make several conclusions. For small /d < 100 pA, a well 
distinguished slope of fog(T) vs. AE corresponding to 
T* w 120 mK is observed when r < 100 s. The typical 
electron temperature in the cryostat Ca is lower, about 
50 mK, so even in this case there are still contributing 
mechanisms beyond the pure TA process. For higher 
barriers and longer hold times, the slope is much lower, 
corresponding to T* 1 K. We presume that the shallow 
slopes are due to PAT at frequency / ~ AE/h, activated 
either by residual photons leaking inside the shields or by 
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FIG. 3. (a) Electrometer data trace at the maximum hystere- 
sis shown in Fig. 1(c) and at low Jd. The trap state changed 
only once at f ~ 14 h. Slow drifts in /d are caused by back- 
ground charge fluctuations of the electrometer, (b) Switching 
trace at the maximum hysteresis with higher I^i > 500 pA. 
(c) Normal state switching trace for sample Sb. 



the back-action of the electrometer. It generates random 
telegraph noise (RTN) whose power spectrum has a rela- 
tively slow decay S{uj) oc 1//^ in the relevant frequency 
scale / » Ici/e (for more details, see Ref. [Ill)- The elec- 
trometer back-action is obviously dominant at the larger 
values of Id, but it cannot be ruled out even in the case 
of low Id and high AE. 

Our main result, very long hold time, was achieved 
with the widest hysteresis loop, i.e. the highest en- 
ergy barrier AE ^ (Aai + Eci) ^ fee x 10 K, and Id < 
100 pA. In that case, time interval between state switch- 
ings was more than 10 hours, which is too long to obtain 
statistically valid data. The example of Fig. 3(a) shows a 
36-hour trace during which the state changed only once, 
14 hours after the start of the measurement. This should 
be compared with the maximum t ~ 20 s for the same 
sample reported in Ref. [llj. At higher Id, see Fig. 3(b), 
the electrometer back-action imposes a significant limita- 
tion for the maximum r. We also did some measurements 
with sample Sb in cryostat Cb • The maximum hold times 
were of the order of 1 h. The difference was most likely 
due to a lower quality of the radiation shielding. 



The trapping times for sample Sb were studied also in 
the normal state in a magnetic field in cryostat Cb- A 
data trace near the maximum of the hysteresis is shown 
in Fig. 3(c). Contrary to the NS hybrid state, where two- 
electron cotunneling is suppressed by the superconduct- 
ing energy gap Q , this leakage mechanism can dom- 
inate in the normal state if PAT is suppressed strongly 
enough. Experimental results are in reasonable agree- 
ment with the expected cotunneling rate corresponding 
to r '--^ 100 s, which we estimated following the approx- 
imations of Ref. [l^ for the case of dissipativc environ- 
ment. For a similar structure without the resistor, the ex- 
pected value for r is of the order of 0.1 ms. We note that 
the cotunneling rates are very sensitive to the parameters 
of the trap, which we have estimated rather indirectly. 
Hence an exact quantitative comparison between theory 
and experiments would require that the parameters were 
determined by measuring the IV curve of the trap, for 
instance, with the help of a cryogenic switch [2|. Never- 
theless, our observations indicate that a high degree of 
shielding of the sample space makes it possible to reach 
and study a quantum-leakagc-limited behaviour of SET 
circuits in the normal state. 

To conclude, a 1000-fold extension of the hold times 
was achieved for a two-junction R-SNS electron trap by 
enclosing it into a specially-designed low-noise environ- 
ment. In the normal state, a switching rate near the 
quantum leakage floor set by the resistively suppressed 
two-electron cotunneling was achieved. The hold times 
of the trap have been measured for the same device over 
a wide dynamic range from 0.1 s to about 10^ s, showing 
high sensitivity to environmental fluctuations, including 
the back-action from the SET electrometer. Device ap- 
plication for sensitive on-chip noise spectrometry seems 
feasible. 
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